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. Abstract. A new method for analyzing the morphological features of 

' point patterns is presented. The method is taken from the study of 

molecular liquids, where it has been introduced for making a statistical 
Q I description of anisotropic distributions. The statistical approach is based 

. on the spherical harmonic expansion of angular correlations. 

> ■ 1. THE METHOD 

o 

In Valdarnini (2001) I propose an alternative method for analyzing the clustering 
morphology which is a generalization of the 2-point function ^(r) (Peebles & 
Hauser 1974) and is based on the spherical harmonic analysis. The method is 
drawn from molecular dynamics simulations, where it has been introduced for 
' studying orientational order of supercooled liquids and metallic glasses (Wang 

& Stroud 1991). Let us consider a system of Np particles. The i — th particle has 
Qh| coordinates r^, in an arbitrary reference frame. For a specified cutoff radius Rc 

Q . all the particles such that |rj — rj \ < Rc are neighbors of i. The line joining i to 

one of the j is termed a bond. The angular coordinates of the vector Aji = fj — fi 
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Qlm{fi)=Y.Ylm{ej,(t>j), (1) 

is the coefficient of the spherical harmonic expansion of the angular density of 
the bonds associated with the particle i. In Eq. |l|, and hereafter, summation is 
understood over all particles j of the distribution such that \ fi — "rjl < Rc- 

The coefficients Qimi^i) are defined as the bond-orientational order param- 
eters and they can be drastically changed by a rotation of the reference systems. 
A natural quantity to consider , which is rotation invariant, is 
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Using the addition theorem for the spherical harmonics, the expression simplifies 
to 
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where Pi is the Legendre polynomial, 

7jfc = cos{0jk) = Aji ■ Aki/{\Aji\\Aki\) 

is the angle between two bonds and the summations in Eq. ^ are then indepen- 
dent of the chosen frame. An useful quantity is the auto-correlation Gi{r) of the 
coefficients Qi{ri)- The function Gi is defined as follows: for all of the Mp pairs 
(i, k), such that |rj — f^l = r it Ar, where Ar is the thickness of the radial bin, 
then Gi{r) is the sum over all of these pairs 

1 A m=l 

This equation can be greatly simplified: let j be the set of neighbors of 
the particle i and p that of the particle k, which satisfy \fj — fi\ < Rc and 
1^ ~ H;! < Rc- Then the summation becomes 

GKO = ^EEEE^Kr,p), (5) 

^^-'p i k j p 

with Tjp being the angle between Aji and Ap^. The summation over the pairs is 
J2i Sfc) with the sum over the particles k only for those particles with \ fi — fk\ = 
r ± Ar. 

The effectiveness of the statistical analysis in quantifying clustering mor- 
phology is studied (Valdarnini 2001) by applying the statistical estimator Gi to 
point distributions produced by an ensemble of cosmological N— body simu- 
lations with a CDM spectrum. The results shown that the statistical method 
defined by the function Gi{r) can be used to analyze the clustering morphology 
produced by gravitational clustering in a quantitative way. The function Gi{r) 
describes anisotropies in the clustering distribution by measuring the degree of 
correlation between the angular densities as seen from two different observers 
separated by r. Gi{r) can then be considered a statistical measure of clustering 
patterns, with different scales probed by varying the input parameters / and 
Rc- With large redshift surveys becoming available in the next few years, the 
proposed statistical method appears as a promising tool for analyzing patterns 
in the galaxy distribution. 
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